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DUDLF 


ABSTRACT 


Fluctuating velocity and temperature data collected from instru- 
ments mounted at three levels on the R/V ACANIA were analyzed. The 
structure of atmospheric turbulence was examined on the basis of 
spectral distributions of wind and temperature fluctuations over the 
ocean. The results of this study were utilized to determine the 
validity of present formulations and prediction techniques for the 
momentum flux and the temperature-structure parameter over water. 

For the stability conditions encountered, the momentum flux 
calculations support present expressions for the constant flux layer. 
Proper directional deviation from prediction of the variation of € 
with height was observed. When corrections for stability were applied, 
proper adjustment toward the predictions were noted. However, com- 
plete adjustment was not evident. This Suggests other factors as 
wind wave effects may be important. C z calculations were in agen 


db 


ment with present formulations and prediction techniques. 
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I. INTRODUCTION 


Interdisciplinary investigations are being conducted at the Naval 
Postgraduate School to correlate optical propagation and atmospheric 
turbulence data in the marine boundary layer. Such information is 
necessary for design of future naval systems. In these investigations, 
personnel from the Meteorology, Mechanical Engineering, Oceanography, 
and Physics departments have been obtaining shipboard observational 
data to evaluate optical propagation in the near-ocean environment. 

The purpose of this study is to examine the structure of atmos- 
pheric turbulence on the basis of spectral distributions of wind and 
temperature fluctuations. Turbulence similarity expressions have been 
tested primarily over land. Only a few over water results have been 
obtained. Observational results from this study will provide a basis 
on which to evaluate and improve existing techniques in predicting 
over water optical propagation. 

During January 1974, an observational experiment was made aboard 
the Naval Postgraduate School research vessel, R/V ACANIA shown in 
Figure 1. This experiment was conducted in Monterey Bay approximately 
two miles offshore. The experimental location provided the near-ocean 
environment for turbulence data collection as well as the platform for 
ship to shore optical propagation data collection. Meteorological 
Measurements of wind and temperature were conducted at three levels. 

The scope of this study concerned the examination of the turbulent 


fluctuations of velocity and temperature in the marine environment. 


lia 
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Figure 1. The R/V ACANIA 
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The examination included an evaluation of the height variation in the 
constant flux layer of (1) € (the viscous dissipation of turbulent 


kinetic energy) and (2) oo (the temperature structure parameter). 
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Il. THEORETICAL CONSIDERATIONS 


A. THE TEMPERATURE-STRUCTURE PARAMETER 
An important parameter in describing optical propagation is the re- 
fractive-index-structure parameter, a - According to Wyngaard, et al 


(1971), c.? 


N can be related to its counter part, the temperature-struc- 


2 : a 
ture parameter, C through a semi-empirical theory. The temperature- 


T 


structure parameter has the form 


Ge. = I@SraOn se) T yr = cay (1) 


a can be estimated by several methods of measurement. Two of these 


methods will be discussed here. 
At a separation distance, r , of the order of inertial subrange 


scales, the temperature structure function in a locally isotropic field 


has the form 


iTGe a Cee = colle (2) 


2 ; : 
C,. can be found by making simultaneous measurements of temperature, 


at two positions separated by the distance r . A second method for 
: : 2 : ; : : 
determination of C.. is by its relation to the one dimensional tem- 


perature spectrum which has the inertial subrange form 


2 gaal3 


i (3) 


o,(K) = BC 


where 8 is a constant and K is the streamwise component of wave 


number. 
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Indirect means of expression C in terms of other meteorological 


T 
parameters is through turbulence similarity theory. Surface-layer- 
similarity theory predicts that when the statistics of the mean and 
turbulent flow fields are properly non-dimensionalized, they become 
universal functions of only a stability parameter. Wyngaard, et al 
(1971) chose to use the Richardson number, Ri, as the stability para- 


meter in doing calculations from overland data. The Richardson 


number is defined by 


_ g/T (d0/dz) 


Ri 
GRICS 


(4) 


where g = acceleration of gravity, T = mean temperature, d@/dz = 
potential temperature gradient, and du/dz = mean wind speed gradient. 
It was shown by Wyngaard that the temperature structure function could 
be expressed as a function of the Richardson number, potential tempera- 


ture gradient and height. 


c. z*/3 (ae/az)* £. (Ri) (5) 


The results of Wyngaard's calculations are shown in Figure 2, verifying 
the similarity theory prediction for unobstructed overland data. 

In order to fully understand the relation of oe to turbulence, 
one must consider € and Ep » the viscous dissipation of turbulent 
kinetic energy and the viscous dissipation of temperature, through 
theoretical relations. Both the balance expression for turbulent 


kinetic energy and Kolmogorov's hypotheses must be examined. These 


theoretical considerations are reviewed in the following paragraphs. 


1S 
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The dimensionless temperature-structure 
parameter versus Richardson Number 


Figure 2, 
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B. BALANCE EXPRESSION FOR TURBULENT KINETIC ENERGY 
The turbulent kinetic energy balance expression has the following 


general form 








2 = 
Vv du ~ Vv 
me i,g-7 - 9 eee 
mame 453 Ox, TOS eee ot ae Wm 7) 
j Oo 3 j 
Zee eT 
-y &, (6) 


Neglecting horizontal eddy transports, energy flux in the vertical can 


be expressed as follows 














The four terms on the right hand side of the equation represent pro- 
cesses which will cause a time variation of turbulent kinetic energy. 
These processes are 


1) Mechanical production of turbulent kinetic energy from the 
mean flow 


ii) Archimedes work term or the energy released due to the effects 
of buoyancy 


iii) Viscous dissipation of turbulent kinetic energy 


iv) Pressure work and flux divergence of turbulent kinetic energy 
due to transfer by turbulent motion 


Further simplification of the kinetic energy balance expression can be 
made by making several definitions and assumptions: 


i) Assuming that the turbulent kinetic energy is in a statisti- 
cally steady state, 
- 2 
i 
7) 


~~ 


a 
ot 





ii) Assuming horizontal homogeneity (i.e., the spatial variations 
of mean fluctuations are negligible), 


-9 -¢ 


Ls 
ox dy 


iii) Specifying the co-ordinate system such that 


iv) Assuming the divergence terms are negligible 


V.V 


O- i P 
OX. 


ns — 
eo 





3 ( 2 


The turbulent kinetic energy balance expression can then be expressed 
in the following simplified form, for the case when the mechanical pro- 
duction of turbulent kinetic is equal to the viscous dissipation plus 


energy released or gained due to buoyance effects 





= ee Son A 
VV x € T v9 (8) 


In cases where the atmosphere is in near=-neutral stability conditions, 
the second term on the right hand side of the equation is negligible 


and Equation (8) can be written as 


2 du 
ue = =€ (9) 
3 


where u,” = -VjV, + The following expression for the vertical gradient 


of the mean wind also applies to the neutral condition 


Se (10) 
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where K is a universal non-dimensional constant (i.e., the Von Karman 
constant) and Z is height. Combining Equations (9) and (10) yields 
an expression to estimate u, based on measurements of e€ at a known 
height Z 


3 


uy, 


€ = oe (11) 


C. VARIATION OF e€ WITH HEIGHT 

Equation (11) can be used to verify the constant flux assumption. 
A straight line plot, with a -l1 slope, of l1n€ against 1nZ implies 
u, is a constant with height. For non-neutral A eae ee in the 
atmosphere, one must retain the buoyancy term in Equation (8). Equa- 


tion (11) can be modified as follows to include the buoyancy effects. 


3 
OK z, Z 
c= xz |  - 7 ag 
or 3 
uy. 
ic [¢., (Ri) ] (13) 
where L = -Tu,?/ek v0 (Monin-Obuhkov Length) 
v9 = heat flux 
Ri = Richardson number 


>) = empirical function of Z/L 
5 = empirical function of Ri 


Equation (13) can be written in the following form 


3 


u 
In € = -In Z + In $,(Ri) + In (—) (14) 
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From overland experiments, various expressions for >, and 5 have 
been derived. For cases of neutral stability, very little variation 
from a -1 slope is noted in a plot of Ine€é versus In Z. In cases 
of instability, an effect on the plot of Ine versus InZ is to 
increase the slope (in a negative sense), to greater than -l. Figure 3 


illustrates these effects. 


D. VARIATION OF sf WITH HEIGHT 


Equation (5) expressed C : as a function of the Richardson number, 


if 


potential temperature gradient and height. The potential temperature 


gradient can be expressed in the following manner 


dé * Z, 
az -z 4a & oo 
where Do v,0/Ku, 


L = Monin-Obuhkov Length 


81 = empirical function of stability 


The Richardson number can also be expressed as a function of Z/L in 


Equation (5) 
Ri=g, ©) (16) 


Combining Equation (5) and Equation (16) results in the following ex- 


2 


pression for C,. 
cee  (o 


7 = Ts a) (17) 
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Figure 3. An Illustration of the Effect of Stability on the 
Variation of e€ with Height 


Zi 





where 3 is an empirical function of stability derived from overland 


data. For cases of neutral stability, 8, = 1 , and 


2-2/3 (18) 


A plot of In C 7 versus ln Z would yield a straight line with a -2/3 


T 
slope ° 
Under stable conditions Wyngaard noted that ee decreased with 
height more slowly than acl - For unstable conditions, overland 


results indicate that 8, approximates the function 


> 4 (2-2/3 


Also, for unstable conditions, Equation (17) can be expressed as 
Zo 4, 2 2/3 ,-4/3 
C,. - 3 oe (Ly Z (20) 


aeplot of In e. versus lin Z would yield a straight line plot with 
a slope of approximately -4/3, assuming overland results are valid 
Over water. Figure 4 shows the predicted variation of cms with 


height for various stability conditions. 


E. KOLMOGOROV'S SPECTRAL SIMILARITY EXPRESSIONS 
The parameters e€ and ex can be estimated on the basis of 
‘Kolmogorov's expressions for universal spectra (in the inertial sub- 


range) of both fluctuating velocity and temperature spectra. The 


inertial subrange in the velocity spectrum represents a range of eddy 
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Figure 4. An Illystration of the Effect of Stability on the Variation 


of Co. with Height 


a3 





sizes where energy is transferred from the production region to the 
Me vination region, (i.e., from lower to higher wave numbers). 
Figure 5 from Lumley and Panofsky (1964) shows the separate regions 
of a typical turbulence spectrum. 

According to Kolmogorov's second hypothesis, the following func- 


tional relations hold 


$(K) = C, (21) 


| 
© 


$7, (K)= 5 & 9 (22) 


where @$(K) and Pp (K) represent one dimensional velocity and tem- 
perature spectra and K represents the wave number. C, and C, are 
empirical constants. er is defined to be eae ; therefore, 
Equation (22) is identical to Equation (3). Equations (21) and (22) 
can be used to find estimates of € and C. by examining respective 
velocity and temperature spectra. 

Turbulent measurements are made at a fixed point in the flow and 
hence spectra are obtained for temporal frequency f . In order to 


use Equations (21) and (22), time and space scales must be related 


using Taylor's (1938) "frozen turbulence" hypothesis 
K = 2nf/u (23) 


where u is the mean wind speed. The hypothesis implies that turbu- 
lence remains unchanged during the time required for it to sweep 


past the probe. The following forms of Equations (21) and (22) feature 
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Figure 5, 


Ek) - - 


Inertial 
Transfer 


Viscous 
Dissipation 


A Schematic Drawing Showing Spectral Transfer and 
Illustrating (A) Energy Containing Region (B) Inertial 
Subrange (C) Dissipation Range and (D) Universal 
Equilibrium Range 
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the transformation from frequency to wave number 


£ $(£) = K $(K) =, 2/3 y-2/3 (24) 
f $,(£)= K $,(K)= C, Gas Ko 2/3 (25) 


A second method for estimating € results from examining the 
upper limit of the inertial subrange of a velocity spectrum. In 
statistical equilibrium, a spectrum of turbulence is independent of 
the source of the turbulence and depends only on € , yY and K 
where y is the kinematic viscosity. The following expression is 


applicable for the inner scale 
3 1/4 
lo = (y?/e)"/ (26) 
where = is proportional to the dissipation wave number. The 
dissipation wave number can be obtained from spectra which show scale 


separation between the inertial subrange and viscous dissipation 


regions. € can then be determined from Equation (26). 
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III. OBSERVATIONS 


A. DATA COLLECTION 

Observational experiments were made from the R/V ACANIA anchored 
in Monterey Bay approximately two miles offshore. Observations were 
made at three levels on the forward mast and bow of the ACANIA. 


Figure 6 shows the mounting arrangements during the experiments. 


B. EQUIPMENT 
1. Fluctuating Velocity and Temperature Measurements 
Velocity sensors were (TSI) Model 1210 "hot wire" probes with 
sliding support shields. Calibration of the sensors was made prior 
to the experiments. Band filtering was applied to the signal with cut- 
off frequencies of 2 hz and 2 khz respectively. The signals were then 
amplified and recorded on a 14 channel Sangamo Model 3562 FM recorder. 
Temperature sensors used were (TSI) Model 1210 platinum wire 
probes with sliding support shields. Electronic conditioning similar 
to that applied to the velocity signals was applied to the temperature 
Signals. The temperature data was analog recorded simultaneously with 
the velocity data. 
2. Mean Velocity, Temperature and Humidity Measurements 
Mean velocity measurements were made at each level using C. W. 
Thornthwaite cup anemometers. The cup revolutions were counted by 
Hewlett Packard 5221A electronic counters. Direct readout of mean wind 
data were then made and recorded in a log. Mean temperature measure- 


ments at each level were made using Hewlett Packard quartz thermometers. 
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Figure 6. 


Mounting Arrangements 
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The signals were passed to a printer, allowing the temperature data to 
be printed out in a time series. Humidity data were collected in a 
fashion analogous to temperature data using a Hydrodynamics Digital II 
system and a printer. 

Figures 7 and 8 show the electronics involved in processing 


both mean and fluctuating data in the ACANIA's laboratory. 


C. CALIBRATION 
1. Calibration of Hot Wire Sensors 
Velocity fluctuation measurements were made with hot wire sensors. 
These velocity sensors were calibrated before each experiment using a 
manometer and TSI 1125 Calibrator in accordance ae procedures outlined 
by Johnston (1974). For any given mean wind speed, there is a corres- 
ponding voltage output. The relationship between wind speed and voltage 


is described by the following equation 
A = a(u)t/2 + b (27) 
where v is voltage and u is the wind speed. Constants a and b 
are calibration curve slope and intercept respectively, and are determined 
experimentally. 
To analyze the recorded velocity data it was necessary to deter~ 


Mine a calibration factor such that 


ce xv' =u! (28) 
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Figure 7, Equipment Used Aboard R/V ACANIA to Obtain Fluctuating 
Velocity and Temperature Data. 
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Figure 8. Equipment Used Aboard R/V ACANIA to Obtain Mean Velocity 
and Temperature Data, 
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Where c = calibration factor, v' 


= voltage fluctuation and u' = 
velocity fluctuation. Differentiation of Equation (27) yields 


avurl? 


a 


0S es (29) 


where c , the calibration is then determined by the relation 


1/2 
Bis AES (30) 


a 


For analysis, determination of the calibration factor consisted of 
(1) obtaining the slope of the calibration curve (2) selecting a mean 
wind speed value observed during the experiment and (3) entering the 
calibration curve with this mean wind value to determine a correspond- 
ing voltage. Figure 9 is an example of an actual hot wire calibration 
curve. 
2. Calibration of Platinum Wire Sensors 

Temperature fluctuation measurements were made using platinum 
wire sensors. A small temperature fluctuation causes a corresponding 
resistance change in the platinum wire. The resistance of a given wire 


can be related to temperature through the following equation 
R = Ro (1 +a (T-To)) (31) 
where R= resistance, T = temperature, To = 0°C, Ro = resistance 


at O0°C, and a = a constant for platinum. The resistance, in ohms, can 


be measured at any given temperature and thus Ro determined for each 


a2 
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wire. Differentiating Equation (31) results in 


dR 


a Ro dT (32) 
or 


R* 


a Ro T! (33) 


which relates a small temperature change to a small change in the 
resistance of the wire. During analysis, resistance fluctuations are 
recorded as voltage fluctuations, through the electronics involved. 

A transfer function, which relates resistance to voltage, must be 


determined during calibration, i.e., 
v' = HR! (34) 


H , the transfer function can be found by substituting a variable 
resistor in place of the wire and recording voltage changes as a 
function of corresponding resistance changes. By plotting these points 
on a graph, the resultant slope will determine the transfer function 

H . Figure 10 is an example of an actual calibration curve used during 
one of the experiments. Substituting Equation (34) into Equation (33) 
results in an expression to relate recorded voltage fluctuations to 


temperature fluctuations 





we (35) 
Or 
T’' = C.v' (36) 


where Cy is the calibration factor. 
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IV. DATA REDUCTION AND ANALYSIS 


A. PRELIMINARY ANALYSIS 

Data period selection was based on evaluations of strip chart 
records of the signals. Criteria for selection of data periods con- 
sisted of (1) signal strength (2) possible averaging time and (3) 
availability of mean profiles. Typical signals meeting the above 


criteria are shown in Figure ll. 


B. DIGITAL METHODS 
1. Analog to Digital Conversion 
Data from selected periods were prepared for further analysis 
by an analog to digital conversion of the magnetic tapes as described 
by Johnston (1974). The analog to digital process was performed on 
the hybrid computer system located on the fifth floor of Spanagel Hall 
at the Naval Postgraduate School. The hybrid system consists of an 
analog computer, COMCOR Ci 5000, interfaced electrically with a 
digital computer, XDS 9300. 
2. Seven Track to Nine Track Conversion 
The analog to digital process yields seven-track octal based 
records written on a digital tape. Analysis of the data was then per- 
formed on - IBM-360/67 computer. This computer requires nine-track 
hexadecimal based digital records. Conversions from seven-track to 
nine-track tapes were made with procedures described in Naval Post- 
graduate School Computer Center Technical Note #0211-08. Two programs 


used for this conversion were TAPEOUT and CONVERT. 
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3. Digital Analysis Programs 


Spectral analysis of the digital tapes was made using FTOR 
(Fast Fourier Transform) and SCOR (Spectral Analysis), programs which 
are available in the IBM 360/67 library. Descriptions of the FFT 
program package have been described by McKendrick (1972) and Johnston 
(1974). The output of th4 SCOR program consists of a tabular computer 
printout as shown in Figure 12 and a graphical plot of the spectrum 


versus frequency as shown in Figure 13. 


C. SCALING PROCEDURES 

During observational and analysis procedures, several gains and 
scaling factors were applied to the recorded signals. Prior to record- 
ing the data, an amplitude gain was applied to the signals. During 
the analog to digital process an additional gain was applied to the 
signal along with a scaling factor. The scaling factor resulted from 
the consideration that the maximum voltage input to the COMCOR Ci 5000 
is 100 volts, and this corresponds to a digital maximum of Dee on the 
XDS 9300. Applying the respective gain and scaling factors resulted 


in a signal input to the IBM 360/67 as follows 


v(volts) x Gl x G2 x SC = v(wolts IBM 360) (37) 
where v = signal from velocity or temperature unit 
Gl = amplification prior to recording 


amplification during A-D process 


@ 
i) 
Il 


SC 


analog to digital scaling factor 
The spectral analysis programs yielded plots of the energy density 


function $(f) , versus frequency f . The corresponding units 
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associated with 9$(f) and f were VOLT“SEC and 1/SEC respectively. 
Multiplying each point of the energy density function by its correspond- 
ing frequency $(f)f results in a quantity with the units of VOLTS“ . 
Taking the square root of this quantity, one can Pare: VOETS ~to 
VOLTS (IBM 360) through Equation (37). Final scaling was accomplished 


by using Equations (29) and (35) which relates volts to either velocity 


Or temperature fluctuations respectively. 
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V. INTERPRETATION OF SPECTRA 


A. SINGLE POINT METHOD 
The Single Point Method is perhaps the most straightforward method 


2 
used to obtain e€ and C,. from velocity and temperature spectra. 


From the spectrum, a wave number, K , within the inertial subrange is 


2 


chosen. Using Equations (21) and (22), calculations of € and C,. 


can then be made directly. 


B. INTERCEPT METHOD 
A second method described by Johnston (1974) can also be used to 


Z ; : 
evaluate c€ and C., from velocity and temperature spectra. This 


method has an advantage over the single point method in that many points 
are used to evaluate the spectral estimate. This is accomplished by 
taking the average slope through many points and determining an inter- 
cept value as discussed in the following paragraphs. 

1. Determination of € 


Equation (24) can be written in the following form 


In KO(K) = -2/3 in K + ln ae (38) 


Plotting 1n K?(K) against -2/3 1n K results in a curve with a slope 


of +1 and an intercept of I1nC a » At the point of intercept, 


iL 


In K equals zero and 


In Kd(K) = In c.e2/ 


1 (39) 
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Letting A equal l1n Kd(K) , the intercept value, € then becomes 


A 
c= )” (40) 
1 


By determining the intercept value, € can be estimated. Figure 14 
shows a typical plot of this type using velocity data at three levels. 
2. Determination of ona 


A similar procedure can be used to determine ee » Equation 


(25) can be written as follows 


In Ko,,(K) = -2/3 In K + In G, oe (41) 


Using the same procedure as described in the determination of € , 


oa can be evaluated at the point of intercept in the following manner 


c 2 - 


T C 2? 


where B= In KO, (K) at the ordinate intercept. Figure 15 is an 


example of an intercept plot using temperature data at two levels. 


C. INNER SCALE METHOD 

A third method of estimating € was to consider the inner scale. 
Equation (26) describes an expression relating € to the separation 
between the inertial subrange and dissipation region of the spectra. 
By choosing the wave number at which the dissipation region occurs in 
the velocity spectra, € can be evaluated. This method was quite sub- 
jective with regard to determining the exact point of scale separation. 
However, the results were qualitatively good when compared to the 


intercept method. 
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VI. RESULTS 


A. SUMMARY OF PERIODS EXAMINED 

Fluctuating velocity and temperature data were collected from three 
levels aboard the R/V ACANIA on 18 January 1974. Selected periods for 
analysis were chosen for three time intervals between 1834 and 1932. 
Table I lists periods examined, heights at which data were collected, 


parameters measured and stability conditions for each period. 


B. VELOCITY RESULTS 
Velocity spectra were obtained for simultaneous measurements at three 


levels during each period. From these spectra € and u, were computed 


* 
by methods described in Sections V-B-1 and V-C. Table II lists ce and 


u calculations for each of these periods. 


x 
The original experimental plans were to measure mean wind data at 

each of the levels. Due to problems encountered with the read-out counter, 

the wind data appeared to be in error. The problem was resolved by use of 

mean wind measurements at one level, the ACANIA's aerovane, and then esti- 

mating a vertical wind field based on the Monin-Obuhkov prediction. A 

Fortran computer program, PROPAR, described by Cardone (1969) was used 

to estimate the vertical wind field. This program uses inputs of (1) 

mean wind speed (2) anemometer height (3) air-sea temperature difference 

and (4) thermometer height to compute the output parameters. These com- 

puted parameters included (1) u, = friction velocity (2) Zo = roughness 


parameter and (3) L = Monin-Obuhkov Length. Mean wind was calculated 


at each level using the following equation 
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TABLE I 
Data Periods Considered in This Study 


18 January 1974 


Time Height Parameter Stability 
(local) (cm) (Ri) 
1834-1854 394.0 Usr,O.Ul TT! -0.136 

(level 1) 
658.0 Uaioe ar 
(level 2) 
822.0 TO), Way el 
(level 3) 
1855-1913 394.0 U,0,0.,U",T’ -0.293 
(level 1) 
658.0 Ueto. ohn” 
(level 2) 
822.0 ne), eae 
(level 3) 
1914-1932 394.0 (ae) (NY oy -0.923 
(level 1) 
658.0 UbsOnUT Ar" 
(level 2) 
$22.0) Usteosue 
(level 3) 
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TABLE II 


€ and Momentum Flux Calculations for 18 January 1974 using 
the Zero Intercept Method 


Height € uy 


2 -3 
(cm sec ~) (en sec l) 


Time 1834-1854 


394.0 83.43 23.60 
658.0 54.74 | 24.33 
822.0 22.83 19.57 


Time 1855-1913 


394.0 342.96 37.80 
658.0 7 2o20 7] rae 
822.0 One 2 22.93 


Time 1914-1932 


394.0 it 2s 26.05 
658.0 13,67 Jee 74 
822.0 7 oa Ap 19.40 
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uU : 

use (1s G) -v, @) (43) 
where Y) is an empirical function of Z/L described by Paulson (1970). 
The purpose in defining the wind field distribution was to obtain advec- 
tion speeds for relating temporal and spatial frequencies. The need for 
wave number (spatial frequency) description was identified in Section 
II-E. 

uineces 16, 17 and 18 show typical velocity spectra at each level for 
one of the periods. In each case a -5/3 slope is observed for the 
inertial subrange region. At higher frequencies, the spectra exhibit 
greater slopes. The regions of different slopes represent different 
scales in the turbulent regime. The scale represented by greater slope 
at higher frequencies, is the dissipation region. Closer investigation 
reveals a decrease in the frequency where the slope changes as height 
increases. However, as mentioned earlier and as can be seen here, the 
exact point of scale separation between the inertial subrange and dissi- 
pation region is difficult to determine. 

Results on € versus height for the velocity data are shown in 
Figures 19, 20 and 21 where Ine is plotted against ln Z for each 
data period. Negative slopes, greater than -l, are noted in each case. 
As stated in Section II-C, this agrees with present boundary layer theory 
for unstable atmospheric conditions. 

The deviation of the slope from -1 was examined to see if itcould 
be accounted for by the influence of the unstable stratification. The 
plots were corrected for stability using Equation (13) where 5 ~ 


(1 - 12Ri) 1/4 - A Slight shift toward a -l slope was observed; however, 
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Figure 16. 
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Figure 17. Velocity Spectrum Level 2 for 1855-1913, 18 Jan 74. 
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Figure 18. Velocity Spectrum Level 3 for 1855-1913, 18 Jan 74. 
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Figure 19. Measured Values of ¢€ vs Height for 1834-1854, 
18 Jan 74, ~ 
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Figure 20. Measured Values of e€ vs. Height for 1855-1913, 18 Jan 74. 
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Figure 21, Measured Values of e€ vs. Height for 1914-1932, 18 Jan 74, 
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complete adjustment to a -l slope was not evident. Similar results 
were found by Johnston (1974). He states "The fact that the stability 
corrections did not completely restore the -1 slope, even in the 
stable case, leads to the examination of the possible effects of wind- 
wave coupling." 

The above results imply that overland prediction for the influence 
of stability are not completely valid over water and that wind-wave 
coupling or wave related velocity fluctuations might be the reason. 
There is still not enough experimental and theoretical knowledge to 
describe the waves’ influence upon the data. 

Several consistent features can be observed in the results. During 
the first two periods there appears to be a bending or "kink" in the 
plots such that they are concave downward. In the third data period, 
the "kink" is concave upward. However, strip chart records show a 
deterioration of the velocity signal at level 2 toward the end of the 
period. This apparently resulted in low estimates of e€ for that 
level. It is not felt that this "kink'’ is due to the ship's influence. 
Welsh (1974) found that ship related motion had little effect on turbu- 
lence measurements at these frequencies. Johnston (1974) did not find 
this characteristic bending in similar plots of lin € versus In Z 
for unstable conditions. In these results, both the intercept method 
and inner scale method exhibit the "kink". Based on only two periods 
of analysis, conclusions as to the cause of this feature cannot be 


asserted. 


C. TEMPERATURE RESULTS 
Fluctuating temperature data were also obtained at three levels for 


each period. At levels 1 and 3, oe was determined from the 
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one-dimensional temperature spectrum as described by Equations (3) and 
(22). At the mid level, two probes were ener. by xr = 10 cm and 

i determined by using Equation (2). Table III lists the C. calcu- 
lations for each of the periods. Temperature spectra at three levels 
are shown for one of the periods in Figures 22, 23 and 24. These are 
typical of the spectra for the other periods examined. 

Results of a versus height are shown in Figures 25, 26 and 27, 
where In ene is plotted against In Z. The general slopes of these 
plots agree reasonably well with the -4/3 prediction for unstable con- 
ditions discussed in Section II-D. Periods 1 and 2 both show close 
agreement with the -4/3 prediction. However, period 3 shows a much 
greater deviation from -4/3. Comparing Richardson numbers from Table I, 
it is noted that period 3 was much more unstable than the first two 
periods. One might expect slightly greater negative slopes for greater 
instability; however, -4/3 should be a limiting slope if overland results 
are valid over water. Since Ce is a function of both € and Eg » 
the possible wind wave effects on e€ discussed in the previous section, 
would also influence cies - This is a possible explanation for any 
deviation in oe from a -4/3 slope. 

Closer examination of ln ee versus 1nZ reveals "kinking" or 
bending similar to that discussed in the velocity results section. In 
both periods 1 and 2, this bending phenomenon is concave downward as it 
‘was for € Pe ations. During period 3 the general appearance of the 
kink in the oy plot approximates the appearance of the kink in the € 
plot, disregarding any influence of a deterioration in the velocity 


Signal during that period. In this case however, the bending is in 


the opposite direction. 


57 





TABLE III 


C : Calculations for 18 Jan 74 using the Zero 
Intercept Method 


Height C 
(cm) eee 
CC scm ) 
Time 1834-1854 
394.0 00312 
658.0 | 00239 
822.0 00168 
Time 1855-1913 
394.0 00600 
658.0 -00558 
822.0 00428 
Time 1914-1932 
394.0 ~ 00934 
658.0 ~ 00416 
822.0 00393 
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Figure 22, Temperature Spectrum Level 1 for 1855-1913, 18 Jan 74. 
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Figure 23. Differential Temperature Spectrum for Fluctuating 
Temperature Data Separated by r = 10 cm, Level 2 
for 1855-1913, 18 Jan 74. 
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Figure 24. Temperature Spectrum Level 3 for 1855-1913, 18 Jan 74. 
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Figure 25. Measured Values of C 5 vs. Height for 1834-1854 
T 
18 Jan 74. 
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Figure 26. Measured Values of C s vs. Height for 1855-1913, 
T 
18 Jan 74. 


63 





- 4.0} 


2 
LnCr 





LnZ 


Figure 27. Measured Values of C,” vs. Height for 1914-1932 
T 
18 Jan 74. 
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The above results give further support to the consideration that 
the bending is actual and not a result of random errors within the 
system, since these spectra were obtained from two independent sources, 
velocity and temperature measurements. However, no further explana- 
tion as to the cause of the bending is offered due to the limited 


amount of data considered. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


Observed height variations of € and cay tend to support present 
formulations and prediction techniques. However, several exceptions 
were noted and will be summarized in the following paragraphs. 

For unstable atmospheric conditions slopes of In € versus In Z 
were greater than -l. When these slopes were corrected for stability, 
there was a tendency for the slopes to shift towards a -l slope. The 
observation was made that complete adjustment to a -1l slope was not 
evident, suggesting other factors such as wind wave effects May be 
important in near-ocean boundary theory. The observation was also made 
that Ine versus In Z was not a linear plot. There appeared to be 
bending or "kinking" which is not predicted in present formulations. 

No definite reason was given to account for the apparent kinking due 
to the limited amount of data examined. 

A plot of ln ch versus ln Z resulted in -4/3 slopes for unstable 
conditions. Some slight deviations were noted and attributed to other 
factors mentioned in the previous paragraph. The characteristic "kinking" 
or bending in the In cos versus In Z_ plots was also evident. 


Spectral analysis appears to be a satisfactory method for determining 


€ , momentum flux and C estimates over the ocean. However, digital 


Te 
methods used are very time consuming due to the amount of tape handling 
and large computer programs used in spectral analysis. A real-time 


spectrum analyzer would increase the amount of data that one could 


examine due to the time saved in the analysis. It is suggested that 
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Subsequent fluctuating velocity and temperature experiments aboard the 


ACANIA be examined using a spectrum analyzer in the analysis process. 
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